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Measuring entanglement is a demanding task that usually requires full tomography of a quantum 
system, involving a number of observables that grows exponentially with the number of parties. 

Recently, it was suggested that adding a single ancillary qubit would allow for the efficient measure¬ 
ment of concurrence, and indeed any entanglement monotone associated to antilinear operations. 

Here, we report on the experimental implementation of such a device—an embedding quantum 
simulator—in photonics, encoding the entangling dynamics of a bipartite system into a tripartite 
one. We show that bipartite concurrence can be efficiently extracted from the measurement of 
merely two observables, instead of fifteen, without full tomographic information. 


Entanglement is arguably the most striking feature 
of quantum mechanics I], defining a threshold between 
the classical and quantum behavior of nature. Yet its 
experimental quantification in a given system remains 
challenging. Several measures of entanglement involve 
unphysical operations, such as antilinear operations, on 
the quantum state Hi, and thus its direct measure¬ 
ment cannot be implemented in the laboratory. As a 
consequence, in general, experimental measurements of 
entanglement have been carried out mostly via the full 
reconstruction of the quantum state JQ. While this 
technique—called quantum state tomography (QST)— 
has been widely used when dealing with relatively low¬ 
dimensional systems 0E1, it is known to become rapidly 
intractable as the system size grows, being outside of ex¬ 
perimental reach in systems with ~ 10 qubits [?]. This 
difficulty lies in having to measure an exponentially- 
growing number of observables, 2 2Ar —1, to reconstruct 
iV-qubits. Such constraint can be relaxed somewhat by 
using, for example, multiple copies of the same quan¬ 
tum state [8], prior state knowledge in noisy dynam¬ 
ics (9), compressed sensing methods 10], or measuring 
phases monotonically dependent on entanglement m- 
However, measuring entanglement in scalable systems re¬ 
mains a challenging task. 

An efficient alternative is to embed the system dynam¬ 
ics into an enlarged Hilbert-space simulator, called em¬ 
bedding quantum simulator (EQS) [12] [13], where un¬ 
physical operations are mapped onto physical transfor¬ 
mations on the simulator. The price to pay, compara¬ 
tively small in larger systems, is the addition of only one 
ancillary qubit and, usually, dealing with more involved 
dynamics. However, measuring the entanglement of the 
simulated system becomes efficient, involving the mea¬ 
surement of a low number of observables in the EQS, in 
contrast to the 2 2N —1 needed with full tomography. 

In this Letter, we experimentally demonstrate an 
embedding quantum simulator, using it to efficiently 
measure two-qubit entanglement. Our EQS uses three 


polarization-encoded qubits in a circuit with two con¬ 
catenated controlled-sign gates. The measurement of 
only 2 observables on the resulting tripartite state 
gives rise to the efficient measurement of bipartite 
concurrence, which would otherwise need 15 observables. 

Protocol. We consider the simulation of two-qubit 
entangling dynamics governed by the Hamiltonian 
H=—go z ®o z , where cr z =|0)(0| — 11)<11 is the z-Pauli ma¬ 
trix written in the computational basis, {|0), 11)}, and 
g is a constant with units of frequency. For simplic¬ 
ity, we let h= 1. Under this Hamiltonian, the concur¬ 
rence [2] of an evolving pure state |i p(t)) is calculated 
as C= \(il>{t)\<J y ® a y K\ip(t))\, where K is the complex 
conjugate operator defined as K\i/j(t))=\ijj(t)*). Notice 
here the explicit dependance of C upon the unphysical 
transformation I\. We now consider the dynamics of the 
initial state |^(0))=(|0) + |l))®(|0)-|-|l))/2. Under these 
conditions one can calculate the resulting concurrence at 
any time t as 

C = \sm(2gt)\. (1) 

The target evolution, e~ lHt \if(0)), can be embedded in 
a 3-qubit simulator. Given the state of interest | if)), the 
transformation 

|^) = |0) ® Re|V>) + |1) <8>Im|^>), (2) 

gives rise to a real-valued 3-qubit state |\f r ) in the cor¬ 
responding embedding quantum simulator. The decod¬ 
ing map is, accordingly, |^>)=(0|4 , )+i(l|W). The phys¬ 
ical unitary gate (7 2 ®ll4 transforms the simulator state 
into a z (gil4|\I/) = |O)(g>Re|'0) — |l)(g)Im|?/>), which after the 
decoding becomes (0|’E)—i(l|\I>)=Re|'0)— *Im|^)=|^*). 
Therefore, the action of the complex conjugate op¬ 
erator K corresponds to the single qubit rotation 
a z ® I4 [T21114] . Now, following the same encoding rules: 
(ip\OK\'i/;)=(Al\(a z —ia x )®0\^), with O an observable in 
the simulation. In the case of 0=(T y ®<j y , we obtain 

C = \{a z ® a y ® <jy) - i(a x ®<Jy®(Jy) |, 


( 3 ) 
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FIG. 1. (a) Qubits 1 and 2 evolve via an entangling Hamil¬ 
tonian H during a time interval f, at which point quantum 
state tomography (QST) is performed via the measurement of 
15 observables to extract the amount of evolving concurrence, 
(b) An efficient alternative corresponds to adding one extra 
ancilla, qubit 0, and having the enlarged system—the embed¬ 
ding quantum simulator (EQS)—evolve via H^ E \ Only two 
observables are now required to reproduce measurements of 
concurrence of the system under simulation. 


which relates the simulated concurrence to the ex¬ 
pectation values of two nonlocal operators in the 
embedding quantum simulator. Regarding the dy¬ 
namics, it can be shown that the Hamiltonian 
H^ that governs the evolution in the simulator is 
<j y ®(ReH)+il 2 <g)(linH) [12]. Accordingly, in 
our case, it will be given by H < ' E ' > =ga y ®<j z ®a z . 

Our initial state under simulation is 
|^(0))=(|0) + |l))G ) (|0)-|-|l))/2, which requires, see 
Eq. § , the initialization of the simulator in 
|’3/(0)) = |0)(g> (|0) + |1)) ( 8 ) (|0) + |1)) /2. Under these con¬ 
ditions, the relevant simulator observables, see Eq. ([3]), 
read {a x ®crcr y ) = sin (2gt) and (a z (g)a y <g)a y )=0, from 
which the concurrence of Eq. (|T|) will be extracted. 
Therefore, our recipe, depicted in Fig. [l] allows the 
encoding and efficient measurement of two-qubit 
concurrence dynamics. 

To construct the described three-qubit simulator dy¬ 
namics, it can be shown (see Supplemental Material) 
that a quantum circuit consisting of 4 controlled-sign 
gates and one local rotation R y ((/))=exp (—ia y 4>), as de¬ 
picted in Fig. [2ja), implements the evolution operator 
U(t)=exp [—ig (u y ®cr z ®cr z ) t], reproducing the desired 
dynamics, with <f> = gt. This quantum circuit can be fur¬ 
ther reduced if we consider only inputs with the ancillary 
qubit in state |0), in which case, only two controlled-sign 
gates reproduce the same evolution, see Fig. [2] (b). This 
reduced subspace of initial states corresponds to simu¬ 
lated input states of only real components. 

Experimental implementation. We encode a three- 












exp |— i 



Ry{4>) 







exp 

|— i 



FIG. 2. Quantum circuit for the embedding quan¬ 
tum simulator. (a) 4 controlled-sign gates and one 
local rotation R y ((j >) implement the evolution operator 

U(t)=ex p , with <j> = gt. (b) A reduced 

circuit employing only two controlled-sign gates reproduces 
the desired three-qubit dynamics for inputs with the ancil¬ 
lary qubit in |0). 


qubit state in the polarization of 3 single-photons. The 
logical basis is encoded according to |/i)=|0), |u)=|l), 
where | h) and |u) denote horizontal and vertical polariza¬ 
tion, respectively. The simulator is initialized in the state 
|’F(0)) = |/i)( 0 )(g) (|/i)W + lu)^ 1 )) ® (| h)^ + |u)( 2 )) /2 of 
qubits 0, 1 and 2, and evolves via the optical circuit in 
Fig. 2 (b). Figure [3] is the physical implementation of 
Fig. 2 (b), where the dimensionless parameter 4>=gt is 
controlled by the angle </>/2 of one half-wave plate. The 
two concatenated controlled-sign gates are implemented 
by probabilistic gates based on two-photon quantum 
interference H5HIZI, see Supplemental Material. 


In order to reconstruct the two three-qubit ob¬ 
servables in Eq. one needs to collect 8 pos¬ 

sible tripartite correlations of the observable eigen¬ 
states. For instance, the observable {cr x ®<jy®<j y ) is 
obtained from measuring the 8 projection combina¬ 
tions of the (|d), |a)}®{|r), |/)}®{|r), |Z)} states, where 
\d)=(\h)+\ v ))/V2, \r)=(\h)+i\v))/V2, and |a) and 1 1) 
are their orthogonal states, respectively. To implement 
these polarization projections, we employed Glan-Taylor 
prisms due to their high extinction ratio. However, only 
their transmission mode is available, which required each 
of the 8 different projection settings separately, extending 
our data-measuring time. The latter can be avoided by 
simultaneously registering both outputs of a projective 
measurement, such as at the two output ports of a polar- 
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FIG. 3. Experimental setup. Three single-photons with wavelength centered at 820 nm are injected via single-mode fibers into 
spatial modes 0, 1 and 2. Gian-Taylor (GT’s) prisms, with transmittance W=1 (i v =0) for horizontal (vertical) polarization, 
and half-wave plates (HWP’s) are employed to initialize the state. Controlled two-qubit operations are performed based on 
two-photon quantum interference at partially polarizing beam-splitters (PPBS’s). Projective measurements are carried out 
with a combination of half-wave plates, quarter-wave plates (QWP’s) and Glan-Taylor prisms. The photons are collected via 
single-mode fibers and detected by avalanche photodiodes (APD’s). 


izing beam splitter, allowing the simultaneous recording 
of all 8 possible projection settings. Thus, an immediate 
reconstruction of each observable is possible. 

Our source of single-photons consists of four-photon 
events collected from the forward and backward pair 
emission in spontaneous parametric down-conversion in a 
6efa-barium borate (BBO) crystal pumped by a 76 MHz 
frequency-doubled mode-locked femtosecond Ti:Sapphire 
laser. One of the four photons is sent directly to an 
avalanche photodiode detector (APD) to act as a trig¬ 
ger, while the other 3 photons are used in the proto¬ 
col. This kind of sources are known to suffer from un¬ 
desired higher-order photon events that are ultimately 
responsible of a non-trivial gate performance degrada¬ 
tion fIM20l , although they can be reduced by decreasing 
the laser pump power. However, given the probabilistic 
nature and low efficiency of down-conversion processes, 
multi-photon experiments are importantly limited by low 
count-rates, see Supplemental Material. Therefore, in¬ 
creasing the simulation performance quality by lowering 
the pump requires much longer integration times to ac¬ 
cumulate meaningful statistics, which ultimately limits 
the number of measured experimental settings. 

As a result of these higher-order noise terms, a sim¬ 
ple model can be considered to account for non-perfect 
input states. The experimental input n-qubit state p exp 
can be regarded as consisting of the ideal state pid with 
certain probability e, and a white-noise contribution with 
probability 1—e, i.e. p exp =£Pid+(l~ s)^ 2 n /^ n - Since the 
simulated concurrence is expressed in terms of tensorial 
products of Pauli matrices, the experimentally simulated 


concurrence becomes C exp =£| sin(2gf)|. 

In Fig. [4j we show our main experimental results from 
our photonic embedding quantum simulator for one cycle 
of concurrence evolution taken at different pump pow¬ 
ers: 60 rnW, 180 mW, and 600 rnW—referred as to 10%, 
30%, and 100% pump, respectively. Figure [4] (a) shows 
theoretical predictions (for ideal pure-state inputs) and 
measured fractions of the different projections involved 
in reconstructing (< 7 x ®a y ®a y ) and (a z ®a y 0<7 y } for 10% 
pump at gt= 7r/4. From measuring these two observables, 
see Eq. we construct the simulated concurrence pro¬ 
duced by our EQS, shown in Fig. [4] (b). We observe a 
good behavior of the simulated concurrence, which pre¬ 
serves the theoretically predicted sinusoidal form. The 
overall attenuation of the curve is in agreement with the 
proposed model of imperfect initial states. Together with 
the unwanted higher-order terms, we attribute the ob¬ 
served degradation to remaining spectral mismatch be¬ 
tween photons created by independent down-conversion 
events and injected to inputs 0 and 2 of Fig. [3] -at which 
outputs 2 nm band-pass filters with similar but not iden¬ 
tical spectra were used. 

We compare our measurement of concurrence via our 
simulator with an explicit measurement from state to¬ 
mography. In the latter we inject one down-converted 
pair into modes 0 and 1 of Fig. [3] For any value of t, 
set by the wave-plate angle <j>, this evolving state has the 
same amount of concurrence as the one from our simula¬ 
tion, they are equivalent in the sense that one is related 
to the other at most by local unitaries, which could be 
seen as incorporated in either the state preparation or 
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FIG. 4. (a) Theoretical predictions (top) and experimentally measured (bottom) fractions involved in reconstructing 

(a x ®<jy®a y ) (left) and (o z (&&y®cry) (right), taken at gt= 7r/4 for a 10% pump, (b) Extracted simulated concurrence within one 
evolution cycle, taken at 10% (blue), 30% (green), and 100% (red) pump powers. Curves represent C=C PP | sin(2gt)|, where C pp 
is the maximum concurrence extracted for a given pump power (pp): Cio%=0.70±0.07, C3o%=0.57±0.03 and Cioo%=0.37±0.02. 
Errors are estimated from propagated Poissonian statistics. The low count-rates of the protocol, see Supplemental Material, 
limit the number of measured experimental settings, hence only one data point could be reconstructed at 10% pump. 


within the tomography settings. 

Figure [5] shows our experimental results for the de¬ 
scribed two-photon protocol. We extracted the con¬ 
currence of the evolving two-qubit state from overcom¬ 
plete measurements in quantum state tomography [4]. 
A maximum concurrence value of C=1 is predicted in 



FIG. 5. Concurrence measured via two-qubit quantum state 
tomography (QST) on the explicit two-photon evolution, 
taken at 10% (blue), 30% (green), and 100% (red) pump pow¬ 
ers. The corresponding curves indicate C=C PP | sin(2grf)|, with 
C PP the maximum extracted concurrence for a given pump 
power (pp): Ci 0 %=0.959 ± 0.002, C 30 %=0.884 ± 0.002, and 
Cioo%=0-694 ± 0.006. Errors are estimated from Monte-Carlo 
simulations of Poissonian counting fluctuations. 

the ideal case of perfect pure-state inputs. Exper¬ 
imentally, we measured maximum values of concur¬ 
rence of Ci O %=0.959 ± 0.002, C 3O %=0.884 ± 0.002 and 
C 10 o% =0.694 ± 0.006, for the three different pump pow¬ 
ers, respectively. For the purpose of comparing this two- 
photon protocol with our embedding quantum simula¬ 


tor, only results for the above mentioned powers are 
shown. However, we performed an additional two-photon 
protocol run at an even lower pump power of 30 mW 
(5% pump), and extracted a maximum concurrence of 
C 5 %=0.979 ± 0.001. A clear and pronounced decline on 
the extracted concurrence at higher powers is also ob¬ 
served in this protocol. However, a condition closer to the 
ideal one is reached. This observed pump power behavior 
and the high amount of measured concurrence suggest 
a high-quality gate performance, and that higher-order 
terms—larger for higher pump powers—are indeed the 
main cause of performance degradation. 

While only mixed states are always involved in ex¬ 
periments, different degrees of mixtures are present in 
the 3- and 2-qubit protocols, resulting in different ex¬ 
tracted concurrence from both methods. An inspection 
of the pump-dependence, see Supplemental Material, re¬ 
veals that both methods decrease similarly with pump 
power and are close to performance saturation at the 10% 
pump level. This indicates that in the limit of low higher- 
order emission our 3-qubit simulator is bounded to the 
observed performance. Temporal overlap between the 3 
photons was carefully matched. Therefore, we attribute 
the remaining discrepancy to spectral mismatch between 
photons originated from independent down-conversion 
events. This disagreement can in principle be reduced 
via error correction mm and entanglement purifica¬ 
tion [23] schemes with linear optics. 

Discussion. We have shown experimentally that entan¬ 
glement measurements in a quantum system can be effi¬ 
ciently done in a higher-dimensional embedding quantum 
simulator. The manipulation of larger Hilbert spaces for 
simplifying the processing of quantum information has 
been previously considered [24] . However, in the present 
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scenario, this advantage in computing concurrence origi¬ 
nates from higher-order quantum correlations, as it is the 
case of the appearance of tripartite entanglement [25ll26j . 

The efficient behavior of embedding quantum simula¬ 
tors resides in reducing an exponentially-growing number 
of observables to only a handful of them for the extrac¬ 
tion of entanglement monotones. We note that in this 
non-scalable photonic platform the addition of one an¬ 
cillary qubit and one entangling gate results in count 
rates orders of magnitude lower as compared to direct 
state tomography on the 2-qubit dynamics. This means 
that in practice absolute integration times favor the di¬ 
rect 2-qubit implementation. However, this introduced 
limitation escapes from the purposes of the embedding 
protocol and instead belongs to the specific technology 
employed in its current state-of-the-art performance. 

This work represents the first proof-of-principle ex¬ 
periment showing the efficient behavior of embedding 
quantum simulators for the processing of quantum in¬ 
formation and extraction of entanglement monotones. 
This validates an architecture-independent paradigm 
that, when implemented in a scalable platform, e.g. 
ions tang, would overcome a major obstacle in the char¬ 
acterization of large quantum systems. The relevance 
of these techniques will thus become patent as quan¬ 
tum simulators grow in size and currently standard ap¬ 
proaches like full tomography become utterly unfeasible. 
We believe that these results pave the way to the efficient 
measurement of entanglement in any quantum platform 
via embedding quantum simulators. 
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SUPPLEMENTAL MATERIAL 

I. Quantum circuit of the embedding quantum 
simulator 

Following the main text, the evolution oper¬ 
ator associated with the embedding Hamiltonian 
= goyf&o,/S>a z can be implemented via 4 control-i? 
gates (CZ), and a single qubit rotation R v {t). These 
gates act as 

CZ ii = |0)(0| (i) 8l (i) + |l)(l| (i) (8) ai j) , (4) 

R l y it) = e~ la y )9t = (cos(gt)I w - i sin (gt)a^ l) ), (5) 

with cr z =|0)(0| —11)(1|, and a y = — i|0)(l|-M|l)(0|. The 
indices i and j indicate on which particle the operators 
act. The circuit for the embedding quantum simulator 
consists of a sequence of gates applied in the following 
order: 


U(t) = CZ 02 CZ 01 R° y {t)CZ 01 CZ 02 . (6) 

Simple algebra shows that this expression can be recast 
as 

U{t) = cos( ff t)I (0) (8) I (1) <8)I (2) - <8 of > 8 f 2) 

= exp 8 of) 8 cr, (7) 

explicitly exhibiting the equivalence between the gate se¬ 
quence and the evolution under the Hamiltonian of in¬ 
terest. 


II. Linear optics implementation 


The evolution of the reduced circuit is given by a R y {t) 
rotation of qubit 0, followed by two consecutive control-Z 
gates on qubits 1 and 2, both controlled on qubit 0, see 
Fig. [6] (a). These logic operations are experimentally im¬ 
plemented by devices that change the polarization of the 
photons, where the qubits are encoded, with transforma¬ 
tions as depicted in Fig. [6] (b). For single qubit rotations, 
we make use of half-wave plates (HWP’s), which shift the 
linear polarization of photons. For the two-qubit gates, 
we make use of two kinds of partially-polarizing beam 
splitters (PPBS’s). PPBS’s of type 1 have transmit- 
tances th =1 and f„=l/3 for horizontal and vertical po¬ 
larizations, respectively. PPBS’s of type 2, on the other 
hand, have transmittances t/,= 1/3 and t v = 1. Their ef¬ 
fect can be expressed in terms of polarization dependant 
input-output relations—with the transmitted mode cor¬ 
responding to the output mode—of the bosonic creation 
operators as 


„t« _ fj~ to 

a p,out — V L P U 


a. 


t(i) 

’p,out 


= 


■p,m + \A - t p a, 

t(0 - yfo 


P®p,in 


t(. i) 

p,in 

tC i) 

p,in ’ 


( 8 ) 

(9) 


where ) stands for the i-th input (output) 

port of a PPBS with transmittance t p for p-polarized 
photons. 

a) 




FIG. 6. (a) Circuit implementing the evolution opera¬ 

tor [7(t)=exp jgcrf ®of )®of ^ , if the initial state 

is |$(0))=|0) (o) ® (|0} (1) + |1> (1) ) ® (|0} (2) + |1} (2) ) /2. 

(b) Dual-rail representation of the circuit implemented 
with linear-optics. Red (blue) lines represent trajectories 
undertaken by the control qubit (target qubits). 

Our circuit is implemented as follows: the first R y {t) 
rotation is implemented via a HWP oriented at an angle 
9 = gt /2 with respect to its optical axis. The rest of 
the target circuit, corresponding to the sequence of two 
control-Z gates, can be expressed in terms of the trans¬ 
formation of the input to output creation operators as 

bhChdh ^ bhChdh 

bhChdy f bfoCfody 

b h c v d h r bfaCydh 

bhCydy t bfiCydy 

b v Chdh t b v c h d h 

byCfi dy i byCfidy) 

by C V d'h ^ byCydh 

by Cy dy ^ b y Cy dy , ( 1 0) 

where b=a^°\ and d=a^ 2 ) denote the creation 

operators acting on qubits 0, 1, and 2, respectively. These 
polarization transformations can be implemented with 
a probability of 1/27 via a 3-fold coincidence detection 
in the circuit depicted in Fig. [6] (b). In this dual-rail 
representation of the circuit, interactions of modes c and 
d with vacuum modes are left implicit. 

The a x and a z single qubit gates in Fig. [6] (b) are im¬ 
plemented by HWP’s with angles 7r/4 and 0, respectively. 
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In terms of bosonic operators, these gates imply the fol¬ 
lowing transformations, 

&x • bh ^ bb v y bfr (11) 

o'z '■ dh —t dh, d v —> d v . (12) 

According to all the input-output relations involved, it 
can be calculated that the optical elements in Fig. ©(b) 
implement the following transformations 

bhChdh -± b h c h dh/ (3v/3) 
b h c h d v -± b h c h d v /{ 3\/3) 
bhCvd h > b h c v d h / (3-\/3) 

bhCydy y bfoCydy / [ 3 \^ 3 ) 

b v c h d h -± b v c h d h /{ 3\/3) 

byCfrdy f byChdy/ (3^) 

byCydh ^ byCydh / 

by Cy dy ^ byCydy / (^ S'vS), (13) 

if events with 0 photons in some of the three output 
lines of the circuit are discarded. Thus, this linear optics 
implementation corresponds to the evolution of interest 
with success probability P = (l/(3\/3)) 2 = 1/27. 

III. Photon count-rates 

Given the probabilistic nature and low efficiency of 
down-conversion processes, multi-photon experiments 
are importantly limited by low count-rates. In our case, 
typical two-photon rates from source are around 150 kHz 
at 100% pump (two-photon rates are approx, linear with 
pump power), which after setup transmission (~80%) 
and 1/9 success probability of one controlled-sign gate, 
are reduced to about 13 kHz (1 kHz) at 100% (10%) 
pump. These count-rates make it possible to run the 
two-photon protocol, described in the main text, at low 
powers in a reasonable amount of time. However, this 
situation is drastically different in the three-photon pro¬ 
tocol, where we start with 500 Hz of 4-fold events from 
the source, in which case after setup transmission, 1/27 
success probability of two gates, and 50% transmission 
in each of two 2 nm filters used for this case, we are left 
with as few as ~100 mHz (~1 mHz) at 100% (10%) pump 
(4-fold events reduce quadratically with pump). Conse¬ 
quently, long integration times are needed to accumulate 
meaningful statistics, imposing a limit in the number of 
measured experimental settings. 


IV. Pump-dependence 

To estimate the effect of power-dependent higher-order 
terms in the performance of our protocols, we inspect the 
pump power dependence of extracted concurrence from 



FIG. 7. Measured concurrence vs pump power. The concur¬ 
rence is extracted from both two-qubit quantum state tomog¬ 
raphy (QST) and the three-qubit embedding quantum simu¬ 
lator (EQS). Straight lines are linear fits to the data. Slopes 
overlapping within error, namely —0.0030 ± 0.0001 from QST 
and —0.0035 ±0.0007 from EQS, show that both methods are 
affected by higher-order terms at the same rate. 


both methods. Fig. [7] shows that the performances of 
both protocols decrease at roughly the same rate with 
increasing pump power, indicating that in both meth¬ 
ods the extracted concurrence at 10% pump is close to 
performance saturation. 

The principal difference between the two methods is 
that in the three-qubit protocol one of the photons origi¬ 
nates from an independent down-conversion event and as 
such will present a slightly different spectral shape due 
to a difficulty in optimizing the phase-matching condi¬ 
tion for both forward and backward directions simulta¬ 
neously. To reduce this spectral mismatch, we used two 
2 nm filters at the output of the two spatial modes where 
interference from independent events occurs, see Fig. [8] 
Note that not identical spectra are observed. This lim¬ 
itation would be avoided with a source that presented 
simultaneous high indistinguishability between all inter¬ 
fering photons. 




Transmission 



A (nm) 


FIG. 8. Spectral filtering. The measured transmission for 
both filters used in our three-qubit protocol qualitatively re¬ 
veals the remaining spectral mismatch. 







